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Abstract —We present the design and measurements of a 90 
GHz prototype of a millimeter-wave channelizing spectrometer 
realized in rectangular waveguide for astronomical instrumenta¬ 
tion. The device was fabricated using conventional high-precision 
metal machining, and the spectrometer can be tiled into a 2D 
array to fill the focal plane of a telescope. Measurements of the 
fabricated five-channel device matched well with electromagnetic 
simulations using HFSS and a cascaded S-matrix approach. 
This motivated the design of a 54-channel i?=200 spectrometer 
that fills the single-moded passband of rectangular waveguide 
in the 130-175 GHz and 190-250 GHz atmospheric windows 
for millimeter-wave spectroscopic mapping and multi-object 
spectroscopy. 

Index Terms —millimeter wave devices, spectroradiometers, 
microwave filters, channel bank filters, spectroscopy 

1. Introduction 

The development and optimization of large format bolo- 
metric arrays for imaging and polarimetry from ground-based 
millimeter-wave and sub-millimeter telescopes has helped to 
revolutionize the fields of cosmology, galaxy evolution and 
star formation. High resolution imaging and spectroscopy of 
individual mm-wave sources is now being done by ALMA 
m, however spectral surveys over wide sky areas and wide 
frequency ranges are not practical with ALMA. The next 
major steps in millimeter-wave imaging and spectroscopy will 
include several science goals. Large area spectral surveys with 
moderate spectral resolution (e.g. R 50 — 200) could be used 
to characterize large scale structure and the star formation 
history of the universe using intensity mapping of emission 
lines such as CO m and CII 0 . Multi-object mm-wave wide¬ 
band spectroscopy with moderate spectral resolution would 
enable galaxy redshift surveys. Further studies of hot gas in 
galaxy clusters through the Sunyaev Zeldovich (SZ) effect 
would be enabled by high angular resolution and moderate 
spectral resolution instruments. 

One of the main new technologies required to achieve these 
science goals is an arrayable wideband spectrometer consisting 
of a compact spectrometer module coupled to highly multi- 
plexable detector arrays. Several groups are working on devel¬ 
oping superconducting on-chip spectrometers (e.g. SuperSpec 
11, El, Micro-Spec (61, DESHIM A (3 and CAMELS (H) 
based on either filter banks or on-chip gratings. Existing mm- 
wave spectrometers in the field include imaging Fourier Trans¬ 
form Spectrometers (91 and the Z-Spec waveguide grating-type 
spectrometer cni. Here we present the design and prototype 
test results for a compact scalable waveguide filter-bank spec¬ 
trometer. This spectrometer can be manufactured with standard 
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high precision machining facilities, and is able to be tested 
both warm and cold independently from the detectors. It can 
be coupled to simple-to-fabricate highly multiplexable kinetic 
inductance detectors (KIDs), or to conventional bolometers. 
The spectrometer is highly complementary to the on-chip 
spectrometers in several ways. First, it could be naturally 
used in the currently undeveloped spectral range from 130- 
250 GHz, which is optimal for measurements of CO line 
emission and the kinetic SZ effect. Also, this device can be 
straightforwardly designed to cover a relatively broad spectral 
resolution compared to the superconducting filters. Finally, 
instead of using KIDs as the detector technology, the device 
could alternately be configured as a room temperature wide¬ 
band backend for mm-wave and cm-wave cryogenic amplifiers 
without the need for downconverting mixers. 

H. Spectrometer Concept and Measurements oe 
THE Prototype 

The design concept of the filter-bank spectrometer is illus¬ 
trated in Fig. which shows a drawing of the five-channel 
prototype filter bank we have fabricated and tested. A direct- 
drilled circular feed horn CD followed by a circular to rect¬ 
angular transition (not shown) couples light from the sky onto 
the main rectangular waveguide. Each channel connects to the 
main waveguide through an E-plane tee that uses an evanescent 
coupling section into a half-wavelength resonating cavity. An 
identical narrow section on the other end of the resonant cavity 
defines the resonating length. The radiation then terminates on 
a detector. The narrow sections of waveguide have a cutoff 
frequency that is 50% higher than the center frequency of the 
passband of the channel. Because the light from the sky is 
below the cutoff frequency of the narrow section, this means 
the section’s impedance is effectively capacitive. However, on- 
resonance the cavity section is effectively inductive, which 
tunes out the capacitive sections, allowing the radiation to 
pass through to that channel’s detector. At frequencies far 
from resonance, this cancellation does not take place, and 
light does not couple to the detector. The center frequency of 
each channel is tuned by adjusting the length of the resonant 
section, and the bandwidth is defined by adjusting the length 
of the narrow capacitive sections. Because the operation of 
the device depends on the narrow coupling sections having 
cutoff frequencies above the maximum frequency incident on 
the device, and because all of our modeling is in the single¬ 
mode limit, the light from the sky will need to be filtered to 
keep the device operating in the single mode limit and prevent 
spurious coupling of high frequencies down the spectrometer 
channels. 
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Fig. 1. Filter concept overview illustrated with the hve-channel WRIO prototype. Light from the feed horn (not shown) comes down the main waveguide, 
and different frequencies are selected off by each of the five channels shown. A device filling the single-moded passband of rectangular waveguide would 
have 54 spectral channels, and has been designed as shown in Fig. The rounded corners are from the machining process, and are treated in the simulation. 
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Fig. 2. Photograph showing the five-channel test device and the measurement 
setup. The source component of a VNA extender was used to generate 
millimeter waves at the input port. The outputs of each spectrometer channel 
were either terminated with a feed horn and absorbing AN-72, or measured 
with a diode detector. This yielded a measurement of the optical efficiency 
and bandwidth of each spectrometer channel. 


To verify this concept, we constructed and tested a prototype 
five-channel spectrometer machined from aluminum for the 
WRIO band. This band was chosen to allow testing using a 
WRIO VNA extender. The dimensions were chosen based on 
simulations in HFSS. Because this spectrometer is realized 
in waveguide, it can operate at both room and cryogenic 
temperatures, which enabled simple and rapid testing. The 
prototype device has one channel with a center frequency fc 
at 80 GHz, three closely spaced channels near 90 GHz, and 
a fifth channel at 105 GHz, each with R = fc/^f ^ 200 
resolution. Just as would be possible with a larger multichannel 
spectrometer, we used E-plane split block construction with 
conventional alignment pins. Nominally no RF current fiows 
across the split. This prototype was machined on a 5-micron 
tolerance CNC milling machine. Tolerances of 1-2 micron 
tolerances that would improve operation at 200-300 GHz are 
regularly achieved on standard high-precision CNC milling 
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Fig. 3. Measured and HFSS-calculated performance of the five-channel 
WRIO test device. The top panel shows the HFSS-simulated passbands of 
the individual channels in the prototype device, shown in colors ranging 
from red to blue. The simulated coupling of the thru detector is shown in 
black. Measured performance is shown in the middle panel, and the table at 
the bottom compares the measured and simulated center frequencies, spectral 
resolutions R, and optical efficiencies (OF). 


machines. 

We measured the passbands of each channel by terminating 
all but one channel using standard gain horns (Quinstar QWH- 
WPRSOO) to dump the power onto absorbing AN-72 foam. 
We then used a broadband square-law diode detector (Pacific 
Millimeter Products WD) to measure the passband of the 
remaining channel, and put a second detector on the thru port. 
We generated millimeter waves using the source component 
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of an OML V10VNA2 VNA extender driven by an Anritsu 
69247A microwave frequency generator. The VNA extender 
produced enough millimeter-wave power to cause the detectors 
to have a non-linear response, so a 20 dB attenuator was 
mounted at the VNA extender source port to keep the detectors 
in their linear response power range. 

The passband of each spectrometer channel was measured 
by sweeping the frequency at the signal generator to generate 
power in 0.025 GHz steps from 75 to 110 GHz, and recording 
the DC voltage on the detector at each step. A detector was 
mounted to each spectrometer channel in succession, while 
terminating all the other channels with the horns. The gain 
variation of the each detector across the passband, as well as 
an absolute power calibration, was established by performing 
a frequency sweep with each detector attached directly to the 
source (with the attenuator still in place). Dividing the spectra 
taken with the detectors attached to the spectrometer by these 
reference spectra gives a measurement of both the absolute 
optical efficiency and the shape of the passband. 

To guide the design of the initial prototype, and for com¬ 
parison with the measurements, we simulated the entire five- 
channel structure in HFSS. The results of the measurements 
and the simulation are shown in Fig. The measured optical 
efficiency of each channel is high, broadly consistent with 
the modeling. The measured center frequencies agree with the 
calculation to better than 0.5%, and the measured spectral 
resolutions agree to within 30%. Out of band coupling is 
observed at the -25 dB level. The standing wave pattern in 
the thru detector suggests an optical path length of roughly 30 
cm, which could suggest a standing wave between imperfect 
terminations in the horns or detectors, and another imperfect 
match somewhere inside the VNA extender. Overall, there is 
good qualitative agreement between the measurement and the 
HFSS simulation, and the selectivities and optical efficiencies 
of the channels are good, suggesting that we can use HFSS to 
design a spectrometer with more channels. 

HI. Cascaded HFSS Simulation Method 

Having verified the spectrometer concept, we designed a 
spectrometer that has a large number of channels spread across 
a wide passband. Here we consider two bands, 130-175 GHz 
and 190-250 GHz. At a spectral resolution of i?=200, 54 
spectrometer channels are needed to fill the single-moded 
passband of rectangular waveguide, assuming we place the 
center frequencies of each channel at the half-power points 
of each of its nearest neighbors. A full 54-channel structure 
would be so large that it would require too much RAM and 
too much CPU time to reasonably simulate in a single HFSS 
run on a workstation. However, we found that cascading the S- 
matrices of single-channel HFSS simulations was as accurate 
as a full simulation, and can be done for a device with a large 
number of channels. We verified the accuracy of this cascade 
method by using it to simulate the five-channel prototype 
device, and comparing that simulation with the full HFSS 
simulation of the entire five-channel structure. 

An overview of the cascade method is shown in Fig. We 
start by drawing the structure of an individual spectrometer 
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Fig. 5. Passbands for the five-channel device simulated with HFSS alone, 
and a cascaded HFSS single-channel simulation for comparison. Since the 
cascade simulation reproduces the details of the HFSS simulation, down to 
-60 dB, this motivates using the method to design a multi-channel device that 
would be too large to simulate with HFSS alone. 


channel, and simulating it in HFSS. The full 3D simulation 
is necessary to calculate the effects of the rounded comers 
left by the machining process, and the fringing fields around 
the waveguide tee. These simulations yield the S-matrix of 
an individual channel. Since the field distribution at each of 
the three ports of an individual channel is nearly identical 
to the fundamental mode of rectangular waveguide, the S- 
matrices of the individual channels can be cascaded together, 
both to each other and to sections of rectangular waveguide, 
to yield an accurate model of the total S-matrix of the 
entire system. We used a function in the scikit-rf module in 
python to perform the cascade, which is an implementation the 
subnetwork growth algorithm ca. The S-matrix of a section 
of rectangular waveguide is composed of simple analytic 
functions, which makes the cascade simulation fast. Cascading 
is much more accurate than naively multiplying the individual 
transmission and refiection values, since it properly treats 
multiple reflections among the internal sub-stmctures in the 
system. 

A comparison of the simulated spectrometer passbands 
using both methods for the five-channel device is shown in 
Fig. We found that we could reproduce all the details of 
the five-channel full simulation down to the -60 dB level by 
using the cascade method. Using this equivalent method for the 
54-channel structure makes it possible to simulate and design 
the device needed for the full instrument. Simulating an entire 
54-channel structure with the cascade method takes only 3 
hours of workstation CPU time to generate the single-channel 
HFSS simulations, and 2 minutes to cascade them together. 

IV. Simulation Results for a 54 Channel Device 

Each channel in the spectrometer has three dimensions: the 
length of the resonator section, the length of the coupling 
sections, and the width of the coupling section. For a 54- 
channel device, that means there are 162 free parameters that 
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Fig. 4. Flowchart illustrating the cascaded S-matrix simulation method for the spectrometer. The blue boxes represent the three-port S-matrices of single 
spectrometer channels simulated using HFSS. The light grey boxes represent the two-port S-matrices of rectangular waveguide that connect the individual 
channels, which are calculated analytically. The S-matrices are cascaded using the scikit-rf module in python. This method successfully simulates the five- 
channel prototype, and is far faster in simulating a multi-channel device than HFSS would be. 


needed to be adjusted to achieve the target center frequency 
and R = 200 resolution for each channel. Tuning all 162 
parameters individually by hand is not feasible, so instead we 
determined the optimal parameters by interpolating between 
hand-tuned designs. Here we illustrate this process for se¬ 
lecting the resonator length dimensions, the coupling section 
lengths and widths were chosen in an identical fashion. First, 
we hand-tuned the dimensions of six channels spread evenly 
across the desired passband by interactively running HFSS 
simulations. The set of center frequencies fdesign of the hand- 
tuned channels was 

/design = { 80*00 GHz, 85.42 GHz, 89.79 GHz,..., 

95.18 GHz, 100.46 GHz, 105.26 GHz} (1) 

and the resonator lengths corresponding to those center fre¬ 
quencies determined by hand-tuning were 

^design = {2.355 mm, 2.068 mm, 1.905 mm, 

1.739 mm, 1.617 mm, 1.523 mm}. (2) 

The coupling section widths were set such that their cutoff 
frequency is 50% higher than the center frequency, and the 
coupling section lengths were adjusted to yield the desired 
R=200 resolution. These sets of hand-tuned dimensions can 
be thought of as a “dataset” of designs with known center 
frequencies and the desired resolution. 

We then scaled these designs to center frequencies across 
the entire band. Putting center frequencies across the whole 
passband and at the half-power points of their nearest neigh¬ 
bors yielded a set of 54 desired center frequencies 

= {80.00 GHz, 80.41 GHz, 80.82 GHz,..., 

103.93 GHz, 104.46 GHz, 105.00 GHz}. (3) 

This can be thought of as Nyquist-sampling the spectral band. 
The hand-tuned designs needed to be scaled somehow to 
determine the resonator lengths to yield these desired center 
frequencies. Since all channels will use the same WRIO 
waveguide, the designs were scaled by the in-guide wave¬ 
length corresponding to each center frequency. For rectangular 
waveguide in the fundamental mode, the in-guide wavelength 
Xg for a frequency / is 



where c is the speed of light, and a is the width of the waveg¬ 
uide. To generate the resonator lengths l}^chn corresponding 
to the desired center frequencies, we linearly interpolated 
between the hand-tuned channels, using the guide wavelength. 
The lengths are therefore 

^NChn design^ ^design") NChn^\ (^) 

where the function mteTp{xdata^ Vdata^ x) linearly interpolates 
between the datapoints {x datadata) to estimate the y value 
corresponding to the input x. The dimensions for the coupling 
section widths and lengths were interpolated similarly from the 
hand-tuned designs. 

Once we had those dimensions, we ran HFSS simulations 
for each single channel. This process was automated using 
the Matlab API for HFSS, and the individual S-matrices for 
each individual channel were stored to disk. This method 
yielded center frequencies that were fairly close to the desired 
values. We took the results of this simulation run as “data” 
and interpolated between them to get designs which are even 
closer to the desired center frequencies. This process, and 
slightly tweaking a few of the channels by hand, yielded a 
set of dimensions where the center frequencies of the final 
54-channel cascaded simulation matched the design goal to 
within an RMS of 0.05% and the spectral resolutions were 
within an RMS of 25% of R = 200. 

To determine the optimal physical spacing between channels 
along the main waveguide, we cascaded the 54 individual 
channel simulations together to yield the 56 port S-matrix that 
simulates the performance of all of the channels in the full 
device. Since changing the spacings only requires repeating 
the cascade step in the simulation, not resimulating in HFSS, 
changing the spacings only takes 2 minutes of workstation 
CPU time to recompute. We were therefore able to simulate 
many spacings and choose the best one. We chose to arrange 
the channels from low to high frequency, and have a constant 
spacing between each channel. Resimulating over a range of 
spacings showed that the optimal channel spacing is 3.065 
mm, which is 3/4 of a wavelength in WRIO guide at 94.2 
GHz. Initial simulations using a different spacing between 
each pair of channels did not yield better performance than 
using constant spacing, but we are investigating this more to 
see if further performance improvements are possible. 

The WRIO design was simulated without any conductor 
loss in the model. The good agreement between the lossless 
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Fig. 6. Simulated passbands of a 54-channel spectrometer. The top panel is 
on a linear scale and shows the 54 Band A passbands which cover 135-170 
GHz and the 54 Band B passbands which cover 190-245 GHz. The black 
curve shows the sum total of all the passbands. The optical efficiency of 
the individual channels ranges from roughly 0.25 to 0.4, which compares 
favorably to the idea-impedance-match case of 0.5. The bottom panel shows 
the simulation down to -60 dB, with selected channels highlighted. Band A 
is indicated on the top x-axis, and Band B is on the bottom x-axis. Out of 
band coupling is simulated to be at the -20 to -30 dB level, or better. 
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Fig. 7. Illustration of a filled focal plane concept, drawn in 3D to scale. Light 
from the telescope comes in from the bottom-right of the illustration, couples 
onto the f/3 feed horns and into the main waveguides, and is selected out into 
frequency channels by the filter bank. KID detectors are in the cards shown 
in green. In this drawing. Band A (the frequencies below the atmospheric line 
at 185 GHz) and Band B (the channels above the 185 GHz line) are fed by 
separate feed horns, but each pixel still only takes up about 10x20 mm of 
focal plane area. We are investigating a diplexer to possibly enable further 
miniaturization, and to keep the system operating in the single-mode limit, 
by feeding both bands with a single smaller horn. 


to Band A gives a total device length of 96 mm, and Band 
B will have a length of 68 mm. These dimensions are small 
enough that fabricating the corresponding detector array cards 
in standard cleanroom processes will be feasible. 


simulations and the measurements of the WRIO prototype 
suggest that in principle we can scale the dimensions and the 
simulation to other slightly higher frequency passbands. Two 
bands lying in atmospheric windows that are interesting for 
observing the SZ effect and CO/CII spectroscopy are the 130- 
175 GHz band, and a 190-250 GHz band. The gap between 
the two bands is to avoid a strong atmospheric absorption 
line. The simulated passbands of the individual channels of 
both a 130-175 GHz band device and a 190-250 GHz device 
are shown in Fig. The passbands of neighboring channels 
cross at the half-power point, which Nyquist samples the 
entire bandwidth. The entire 130 GHz to 250 GHz passband 
does not fit in the single-moded bandwidth of a rectangular 
waveguide, so we will use either a single horn with a diplexer, 
or two independent feed horns, to place the full bandwidth 
into a lower Band A and an upper Band B, on either side 
of the atmospheric line at 180-185 GHz. Since the device is 
small enough to fit under the footprint of the feed horn, a 
linear array of these spatial pixels is formed in one direction, 
which all feed a single wafer of KID detectors. These linear 
arrays can then be tiled in the other direction, like vertical 
cards in a motherboard, to form a filled focal plane array of 
spectrometers. 

A concept schematic of a 2 x 2 array of spatial pixels is 
shown in Fig. [ 7 ] For a single spatial pixel in the focal plane, the 
light will come in from the horn, down the main waveguide, 
and be directed into the 54 individual spectral channels. In the 
prototype WRIO design, the physical center-to-center spacing 
between channels is a constant 3.065 mm. Scaling the design 


A. Loss and Machining Tolerances 

Scaling the lossless simulation to design higher frequency 
devices is a valid approximation as long as conductor loss 
continues to be negligible in the higher passbands. For a rect¬ 
angular waveguide, the attenuation constant due to conductor 
loss is 


OLn 




b “1“ 2ft' 






( 6 ) 


where a and b are the width and height of the waveguide, a 
is the metal conductivity, and r] is the free space impedance 
ca. Since the resonator quality factor is defined in terms of 
fractional energy lost per cycle, the Q due to conductor loss 
is approximately 


Qi 
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( 7 ) 


This loss Q will degrade the actual spectral resolution Rtot 
from its nominal lossless design value R by 


Rtot R Qioss 

For the WRIO prototype with aluminum’s room temperature 
DC conductivity (141, the calculated attenuation constant of 
aluminum WRIO waveguide at 105 GHz is 0.25 m“^. This 
implies a calculated Qioss of approximately 7,000, which is 
high enough to explain the good performance of our prototype. 
Since our prototype was roughly R = 200, directly measuring 
the impact of this loss Q would have required fabricating and 
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Fig. 8. Calculation of the limiting spectral resolution due to conductor loss 
and machining tolerance. The top three curves show loss Q calculated with 
literature room temperature conductivity values for OFHC copper, gold, and 
aluminum, which range from roughly 2,000 to 8,000. The bottom curve shows 
the highest resolution for which nearby channels could be distinguished with 
a 1 micron tolerance milling machine. This shows that a i?=200 spectrometer 
would not be limited by either effect at operating frequencies below roughly 
700 GHz. 


testing another test device with high design R, which we have 
not yet done. 

Scaling up the maximum operating frequency, and the corre¬ 
sponding reduction in the waveguide dimensions, both increase 
the loss. However, for a maximum operating frequency of 250 
GHz, the loss Q is calculated to be roughly 4,500, which is still 
high enough that it would not be expected to limit the perfor¬ 
mance of a = 200 device. Sputter coating a superconducting 
niobium layer onto the device would eliminate conductor loss 
below ^ 700 GHz. 

Another factor that will limit the performance of the 
spectrometer is machining tolerance. If the dimensions of 
nearby channels in the spectrometer differ by less than the 
mechanical tolerance to which they are fabricated, then those 
channels cannot be resolved. This means that the maximum 
practical resolution of our half-wave resonator channels is 
approximately 

TD _ __ /Q\ 

Machine Tolerance' 

At ASU, and at other high-precision machine shops, we have 
a precision milling machine that regularly achieves a 1 micron 
tolerance, which means that as shown in Fig. a i? = 200 
device should not be limited by machining tolerance below 
about 700 GHz. 

V. Conclusions 

Compact spectrometers that can be tiled into focal plane 
arrays are an important enabling technology for the next 
generation of millimeter-wave and sub-millimeter astronomy. 
Complimentary to the other technologies currently under de¬ 
velopment, the measurements and modeling of the waveguide 
filter-bank spectrometer presented here show that it is a 


promising approach for future instruments. Measurements of 
the prototype show that cascading the S-matrices from HFSS 
simulations of individual channels is a good way to model 
this class of device. This enabled us to design a full R=200 
spectrometer that fills single-moded passband of rectangular 
waveguide. We are currently scaling up our prototype to a five- 
channel test device for the WR5 band (160-210 GHz). We will 
then test it with another VNA extender at room temperature 
to verify that conductor loss and machining tolerances do not 
limit the performance at higher frequencies. We then plan 
to fabricate and test a full 54-channel device with cryogenic 
detectors for the 130-175 GHz and 190-250 GHz bands to 
verify that the optical efficiency and spectral performance are 
as good as the modeling predicts. 
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